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Abstract. I present a new method to unveil the history of cosmic accretion and the 
build-up of Supermassive Black Holes (SMBH) in the nuclei of galaxies, based on ob- 
servations of the evolving radio and (hard) X-ray luminosity functions of AGN. The 
fundamental plane of black hole activity discovered by Merloni, Heinz & Di Matteo 
(2003) is used as a mass and accretion rate estimator. I adopt the local BH mass func- 
tion as a boundary condition to integrate backwards in time the continuity equation for 
the SMBH evolution, neglecting the role of mergers. Under the most general assump- 
tion that accretion proceeds in a radiatively efficient way above a certain rate, and in a 
radiatively inefficient way below, the redshift evolution of the mass and accretion rate 
functions are calculated self-consistently. The only tunable parameters are the accre- 
tion efficiency and the critical ratio of the X-ray to Eddington luminosity at which the 
transition between accretion modes takes place. The evolution of the BH mass function 
between z — and z ~ 3 shows clear signs of an anti-hierarchical behaviour: while the 
majority of the most massive objects {M > 10^) were already in place at 2 ~ 3, lower 
mass ones mainly grew at progressively lower redshift. As an example, I will discuss 
the consequences of these results for the lifetimes of accreting black holes. 

1 Introduction 

It has been known for the last ten years that the cosmological evolution of 
massive galaxies shows signs of 'down-sizing' i.e. of a progressive decrease 
of the typical mass of actively star forming galaxies. Many pieces of evidence, 
brought forward also during this meeting (see e.g. the contributions of Bender, 
Kauffmann and Danese to these Proceedings), suggest that the baryonic physics 
of star formation processes counterbalance the hierarchical (bottom-up) growth 
of dark matter structures, favouring the early formation and growth of the most 
massive (and dense) stellar aggregates. 

The ubiquity of SMBH in galactic nuclei, and their tight relation with their 
hosts' bulge properties |2j seem to indicate that the formation and growth of 
galaxies and of their resident black holes have proceeded in parallel, and probably 
influenced each other in the early stages of galaxy formation. As a matter of 
fact, the number of theoretical studies dedicated to AGN feedback in galaxy 
formation has flourished in the last five years (see e.g. Loeb's contribution in 
this proceedings, and references therein). Furthermore, a recent comprehensive 
study of 23 000 local AGN carried out by the Sloan Digital Sky Survey (SDSS, 

have demonstrated, in a direct quantitative way, that accretion onto SMBH 
and formation of stars are tightly coupled even in the local universe. 
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Is it possible to learn more about the formation and growth of structures 
by just looking at the evolution of the AGN population? The aim of the work 
presented here is to show to what extent this is indeed possible, and to describe 
a robust, self-consistent way to unveil the history of cosmic accretion and the 
build-up of SMBH in the nuclei of galaxies, in the form of their evolving mass 
function. The methodology and the main results will be discussed in the next 
section, while in section |3| I will expand on the consequences of these results for 
the issue of QSO lifetimes. Section 4 will summarize my conclusions. 

2 The Evolution of the SMBH Mass Function 

Under the standard assumption that black holes grow mainly by accretion 0], 
the cosmic evolution of the SMBH accretion rate and its associated mass density 
can be calculated from the luminosity function of AGN: (j){Lhoi,z) = dA^/dLboi, 
where Lboi = ^mdMc^ is the bolometric luminosity produced by a SMBH ac- 
creting at a rate of M with a radiative efficiency Crad- In practice, the accreting 
black hole population is always selected through observations in specific wave- 
bands. Crucial is therefore the knowledge of two factors: the completeness of any 
specific AGN survey, and the bolometric correction needed in order to estimate 
Lboi from the observed luminosity in any specific band. On both these issues, 
huge progress has been made in the last few years (see e.g. [H])- 

In order to progress from the study of BH mass densities to that of BH 
mass functions, we need to break the degeneracy between mass and accre- 
tion rate of any given observed AGN. While in most semi-analytic works this 
is done by assuming a constant Eddington ratio for all sources, here I will 
propose an alternative, physically motivated, method. In a recent paper |H] 
it has been shown that the hard (2-10 keV) X-ray luminosity of an accret- 
ing black holes is related to its mass and its core radio (at 5GHz) luminos- 
ity by the following relation (the "fundamental plane" of black hole activity): 
logiR = (0.60tg:")logix + (0.78tg:Jg)logM + 7.33^4^^; which can be in- 
verted to relate BH masses to observed nuclear radio and X-ray luminosities: 
logM ~ (;(logLR, logLx)- One of the consequences of this relation is that, in an 
ideal case, the conditional radio/X-ray luminosity function of active black holes, 
i.e. the number of sources per unit co-moving volume per unit logarithm of radio 
and X-ray luminosity, \1^c{Lr, Lx), could be used to reconstruct the mass func- 
tion of the underlying black hole population. In fact, the current lack of the exact 
knowledge of !f'c(-^R,^x) can be (at least partially) superseded, given the two 
separate radio, (/)fi(Lr, z), and X-ray, (pxiLx, z), luminosity functions at redshift 
z, and an independent estimate of the black hole mass function, (j)yi{M, z) at the 
same redshift. By taking into account the fundamental plane relationship, we 
have that, at any z, the conditional luminosity function •Z'c has to satisfy the 
following integral constraints: 

/•oo 

(/.,(L,)dlogL, = / >Fc{U,Lj)d\ogLj (i,j) = (R,X) (1) 
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(t>MiM)d\ogM 



11 ^ciLx,LB.)d\ogLnd\ogLx. (2) 

J J log M <g<log M+d log M 

Given observational estimates of (f>x, fpR ^ind (f>M, we start with an initial guess 
for tf'ci and proceed via successive iterations, minimizing the differences between 
the projections of the conditional luminosity function onto the X-ray and radio 
luminosity axes and the observed luminosity functions, until a conditional LF is 
obtained simultaneously satisfying eqs. and Once such an estimate of tf'c 
is found, it is possible to derive the local distribution of the X-ray to Eddington 
ratio, and from this, given an appropriate bolometric correction (i.e. a specific 
functional form Lx ~ Lx{M, m), where m = eMc^ / L^^d), the desired accretion 
rate function. 

The redshift evolution of the SMBH population can then be computed inte- 
grating backwards the continuity equation that describes SMBH evolution driven 
by accretion only [Jj: 

dMM,t) d[MM,t)-{M{M,t))] _ 

dt + dM 

where the mean accretion rate as a function of black hole mass and time, (M) 
is calculated directly from the accretion rate distribution function at time t. 
Starting from z = 0, the local BHMF, as determined independently from the 
galaxy velocity dispersion distribution and the M — a relation, can be used 
as a boundary condition to integrate cq. ||2Jl up to the redshift where hard X- 
rays and radio luminosity functions of AGN can be reliably estimated. The 
only parameters needed are the accretion efficiency e, and the functional form 
Lx^LxiM, m). 

Here I present the results obtained using the most recent determinations of 
the hard X-rays and radio LF |9llUj . fixing e = 0.1, and assuming for the function 
Lxifn-) a simple broken power-law expression, with Lx oc fn? at low accretion 
rates (radiatively inefficient regime) and Lx oc to"-^ at high luminosities (radia- 
tively efficient regime, see |5I8| for details), with the transition being placed at 

= 10-3 HI]. 

The evolution of the black hole mass function is shown in the left panel of 
Figure ^ As opposed to the standard picture of hierarchical mass build up of 
dark matter halos in GDM cosmologies, supermassive black holes growing by 
accretion between z 3 and now have a mass function which is more and more 
dominated by larger mass objects at higher redshift. Thus, most of the more 
massive black holes (M > 10^) were already in place at z ^ 3, and as a result, 
the "typical" SMBH mass decreases with decreasing redshift. SMBH appear to 
be growing in anti- hierarchical fashion 5 12. 

The right panel of Figure^shows instead the redshift evolution of the derived 
accretion rate function (expressed here as X-ray to Eddington ratio, ix/-^Edd)- 
While the number of sources accreting at low rates increases monotonically with 
decreasing redshift, the situation is different for rapidly accreting objects, that 
increase in number with increasing redshift. The cut-off redshift, above which the 
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Fig. 1. Left Panel: Redshift evolution of the SMBH mass function (BHMF), from 
redshift 3 (lower curve) till redshift 0.1 (upper curve). Right Panel: Redshift evolution 
of the SMBH accretion rate function (expressed here in terms of X-ray to Eddington 
luminosity ratio), from redshift 3 till redshift 0.1; The vertical dashed line marks the 
adopted value of the critical accretion rate where a transition occurs between radiatively 
inefficient (below) and efficient (above) accretion. 



number of sources declines again, is a function of the typical X-ray to Eddington 
ratio, being lower for lower accretion rate sources. Once again, it is clear that 
the "typical" X-ray to Eddington ratio, which can be approximately identified 
with the knee of the accretion rate functions of Fig.l, decreases with decreas- 
ing redshift, crossing the critical rate that separates radiatively inefficient from 
efficient regimes (here fixed at 10~^, vertical dashed line) around z ~ 0.5. 

3 Lifetimes of Active Black Holes 

In all models that derive the properties of the SMBH population from the ob- 
served QSO evolution, a key element is represented by the typical quasar lifetime 
or by the almost equivalent activity duty cycle. However, the significance of these 
parameters is limited to the standard case in which, on the basis of an observed 
luminosity function in a specific waveband, one tries to derive the distribution 
of either BH masses or accretion rates. Usually, a constant Eddington ratio is 
assumed in this case, which implies that QSO are considered as on-ofF switches. 
Then, the duty cycle is simply the fraction of black holes active at any time, and 
the lifetime is the integral of the duty cycle over the age of the universe (see e.g. 
[T^ and references therein). 

The picture discussed here is different, in that a broad distribution of Edding- 
ton rates is not only allowed, but actually calculated for the SMBH population at 
every redshift. When this is the case, a more meaningful definition of activity life- 
time is needed. Let us first define the mean Eddington rate for object of mass Mq 
at redshift z = (m(Mo,2:)) and then introducing the mean accretion weighted 
lifetime of a SMBH with a given mass today: t{Mo,z) — J^{m{Mo, z'))jp-dz' . 
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Fig. 2. Partial mean accretion weighted lifetimes (in years) of SMBH with mass today 
Mo, calculated for three different redshift intervals: < z < 1 (solid line), 1 < z < 2 
(dashed line) and 2 < z < 3 (dot-dashed line). The horizontal dotted line is the Salpeter 
time for accretion efficiency of 10%. 



The ratio of t(Mo, z) to the Salpeter time, is = eMc'^/L^dd = (e/0.1)4.5 x 10^ 
yrs, gives the mean number of e-folds of mass growth for objects with mass Mg 
up to redshift z. The ratio of t{Mo,z) to the Hubble time iHubbic(-z) ~ H{z)^^, 
instead, is a measure of the activity duty cycle of SMBH. It is also interesting 
here to calculate "partial" lifetimes in a given redshift interval Az — (zi, zj): 

dt 

At{Mo,Az) = J {m{M„,z')) — dz'. (4) 

In Figure [3 I show At{Mo, Az) for three redshift intervals: < z < 1; 
1 < z < 2 and 2 < z < 3. The accretion weighted lifetime for BH of any 
given mass between < z < 3 is of course just the sum of the three. The 
anti-hierarchical nature of mass build-up in actively accreting AGN and QSOs 
is again clearly illustrated by this plot. In fact, the major growth episode of a 
SMBH must coincide with the period when At > ts- This happens at z < 1 for 
Mo < 10'^•^ between redshift 1 and 2 for lO^ '^ < Mq < 10^ '^, and at 2 < z < 3 
for 10^'^ < A/q < 10^'^. Supermassive black holes with masses larger than Mq ~ 
10^'^ today, must have experienced their major episodes of growth at redshift 
higher than 3. It also interesting to note that the objects that dominate the 
SMBH mass function today, i.e. those in the range of masses around lO^'^M©, 
where Mo(/)m(-^Oj z — 0) peaks, mainly grew around z ~ 1, which is when most 
of the X-ray background light we see today was emitted . 
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4 Conclusion 

I have presented a new method to study the growth of accreting supermassive 
black holes, based on the simultaneous evolution of the AGN radio and hard (2- 
10 keV) X-ray luminosity functions. The method is based on the locally observed 
trivariate correlation between black hole mass, X-ray and radio luminosity (the 
so-called fundamental plane of black hole activity, [HI)- Thanks to this correla- 
tion, it is possible for the first time to break the degeneracy between luminosity, 
mass and accretion rate: QSO and AGN not only grow in mass during their 
evolution, but also accrete a different rates depending on their mass and age. 

Qualitatively (i.e. independently on the values of the model parameters), 
the evolution of the black hole mass function between z = and z ~ 3 shows 
clear signs of an anti-hierarchical behaviour: while the majority of the most 
massive objects (M > 10^) were already in place at z 3, lower mass ones 
mainly grew at progressively lower redshift, so that the average black hole mass 
increases with increasing redshift. At the same time, the typical accretion rate 
of SMBH decreases with decreasing redshift. Consequently, the lifetime of an 
actively growing BH, measured through the partial mean accretion weighted 
lifetime (see sectionlS)), is a strong function of both redshift and black hole mass. 

Broadly speaking, what is presented here is the SMBH analog of the down- 
sizing of star forming galaxies. Future more direct and detailed comparisons be- 
tween SMBH mass functions, such as those shown here, and galaxy luminosity 
functions, both at low and high redshifts, will hold the key of our understand- 
ing of the complex interplay between formation and growth of black holes and 
galaxies. 
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